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Definition 
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m2 
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Coupling function 

dimensionles s 

f 

Frequency of fluctuation 

Hz 
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Constant complex number 
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Normalized pressure fluctuation 

dimensionles s 
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rev./ s, (rps) 
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dimens ionle s s 
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Pressure 

N/m2 

Ap 

Deviation of pressure from time-mean value 

N/m2 
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Flow rate through turbopump 

m^ / s 

AQ 

Amplitude of flow fluctuation 
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Radius from axis of rotation 
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Constant for conformal transformation 
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Spacing of airfoils in cascade 
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S(co, 0-) 

Sinusoidal gust function 
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U, V 

Induced velocity 

m/ s 
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Moving velocity of rotor cascade 

m/ s 
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Through-flow velocity parallel to airfoil 
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Absolute velocity 

m/ s 

W 

Unsteady inlet flow in cascade direction 

m/ s 
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Relative velocity to rotor cascade 

m/ s 
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Coordinate, parallel and perpendicular 
to airfoil 
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X, Y 

Coordinate, axial and cascade direction 

m 

Y 

Distributed circulation 

m/ s 

r 

Cir culation 

m2 /s 
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Definition 

Units 

0 

Polar angle 

radian 

X. 

Stagger angle of cascade 

radian 
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Angular velocity = 2 tt f 

radian/ s 
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Density of working liquid 

kg/m^ 

cr 

Solidity of cascade = f / s 
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‘I’ 

Flow coefficient 
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Velocity potential 

m^ /s 


Pressure rise coefficient 

dimensionless 
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Reduced frequency 
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Gravitational potential 

m^ /s^ 
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Pertinent to translatory oscillation 
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Pertinent to sinusoidal gust oscillation 
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SUBSCRIPTS 

Time- mean value 


1 

Upstream of cascade 


2 

Downstream of cascade 


a 

Axial component 


C 

Due to conduit effect 


d 

At pump delivery port 


qs 

Quasi- static condition 


R 

Rotor cascade 


s 

At pump suction port 


S 

Stator cascade 


R + S 

Due to pumping action of rotor and stator 


u 

Component parallel to cascade direction 


W 

Pertinent to wake 
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ANALYTICAL AND EXPERIMENTAL STUDY OF 
DYNAMIC CHARACTERISTICS OF TURBOPUMPS 


SUMMARY 

The response of pressure rise of turbopiimps to fluctuating flow rate, 
termed the ’’dynamic characteristics, ” was studied analytically and experi- 
mentally. Unsteady flow around a two-dimensional, linear cascade of air- 
foils was solved and the frequency response of the deflection angle of the 
cascade was calculated for the periodically oscillating inlet condition. Based 
on the imsteady cascade theory, the dynamic characteristics of turbopumps 
were analyzed. As a sample of the theoretical application, an exact method 
for calcTolating the dynamic characteristics of a single-stage axial flow pump 
is shown. For the evalxoation of an approximate response, a simplified model 
theory was also developed, which enables the calculation to be made without 
entering into the details of the pirmp design. 

An experimental investigation was made to determine the dynamic charac- 
teristics of a centrifugal pump. The test showed that the analytical calculation 
was close to the experimental result. The limiting fluctuating frequency under 
which the dynamic characteristics agree with the quasi- static ones, was deter- 
mined from the test result. The time constant of the response of pressure 
rise to a step-like change of flow rate was also derived. 

The entire study was made imder incompressible and non- cavitating con- 
ditions. If the cavitation occur s in a turbopump, the phenomena change drasti- 
cally from those described in this study. 

INTRODUCTION 


The pressTire rise through a turbopump^i^ varies with the flow rate. The 
relationship between them is indicated by the characteristics curve of the 
pump. This characteristics curve is valid only when the pump operates in 
the steady- state condition, that is, with a constant flow rate. When the flow 
rate is no longer steady and changes rapidly, the pressure rise cannot 
respond quickly enough to follow along the steady- state characteristics curve. 
This results in a considerable change in the actiial characteristics c\u:ve. The 
purpose of this study was to determine analytically and experimentally the 
correlation between the flow rate oscillation and the resulting pressure rise 


* Turbopumps include axial flow pumps, mixed flow pumps, centrifugal 
pumps etc. , which are all classified as turbomachinery. 


oscillation of tiirbopumps, termed dynamic characteristics. 

Since the unsteady condition is restricted (in this study) to the case in 
which the flow rate oscillates sinusoidally around its mean value with a 
constant frequency, the dynamic characteristics are identical to the frequency 
response of pressure rise to the flow rate. The dynamic characteristics are 
especially important for analyzing the stability of a system in which the 
pump is installed. Figure 1 shows a tiirbopiomp whose flow rate, Q, oscillates 
periodically around its mean value, Qq, because of some outer disturbances. 
The pressure at the suction port, pg, and at the delivery port, p^j, fluctuate also 
according to the fluctuating flow rate, and the pressure rise through the pump, 
Pd " Ps» fluctuates correspondingly. The amplitude and phase correlations 
between the pressure rise, pj - Pg, and the flow rate, Q, were the essential 
information pursued in this study. The rotational speed of the pump was 
assumed to be constant and unaffected by the fluctuating operation. 

Many studies (Reference 1 to 3 )have been made to predict analytically, the 
steady state characteristics curves of turbopumps. Most of these are derived 
from the consideration of the flow field aroundthe cascade of airfoils, because the 
behavior of the cascade is the principal determinant of the characteristics of' 
the turbopump. For the same reason, the behavior of the cascade is the 
determining factor of the dynamic characteristics of a turbopump. There- 
fore, an analysis of the flow through a cascade of airfoils Tinder the oscillating 
inlet flow condition was made as the first step in the theoretical approach and 
is described in Section I of this paper. Using this basic result, another analy- 
sis was performed to determine the dynamic characteristics of turbopumps 
and is described in Section II. Experimental research was also conducted to 
determine the dynamic characteristics of a centrifugal pump. The test results 
and the comparison with the analytical results are described in Section III. 
Interest in this study is limited to the incompressible, noncavitating case. 
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SECTION I. ANALYSIS OF OSCILLATING FLOW THROUGH A CASCADE 


Two-dimensional, incompressible, inviscid, and noncavitating flow 
around a cascade of airfoils was analytically determined for a periocically 
oscillating inlet flow. Figure 2(a) illustrates the inlet condition to the cascade. 
The unsteady inlet flow, Vj, can be considered as a superposition of a steady 
inlet flow, Vjo> an unsteady flow, w(X, Y, t). The unsteady flow, w, can 

be further divided into two unsteady components, w* and w**, parallel to the 
direction of cascade and to the direction of zero-lift, respectively. 

The unsteady component, w**, causes no change of circulation around the 
airfoils, and represents an oscillating translation flow that can be simply 
superimposed on the flow solution. On the other hand, the unsteady compo- 
nent, w*, parallel to the cascade direction produces a related imsteady circu- 
lation around the airfoils, which is to be determined in this Section. From 
the previous consideration, the analysis was conducted specifically for the 
case in which the unsteady component of the inlet flow is parallel to the cas- 
cade direction, as illustrated in Figure 2(b), without losing any generality. 

The unsteady component, w, now parallel to the cascade direction, is generally 
a fvinction of time and the coordinates. 

There are fundamentally two types of unsteady components in the inlet 
flow. One is the case in which the unsteady component, w, is expressed as 
w = w Exp (ivt)*. In this case, termed “ translatory oscillation,” w is instan- 
taneously uniform in the region preceding the cascade and oscillates periodi- 
cally with the frequency f = v/27t. Translatory oscillation is associated with 
the oscillation of the flow rate, that is, oscillation of axial velocity, In 

the other case, w is expressed as w = w Exp (ivt - ivX/Vg^)*. This oscilla- 
tion represents a sinusoidal gust propagating downstream with a velocity of 
Vg^. This case is termed '‘sinusoidal gust oscillation” and occurs as the result 
of unsteady circulations of airfoils in the preceding cascade. 

To avoid excessive mathematical complexity, the unsteady flow was 
analyzed for cascades of flat plates having a zero stagger angle for two types 
of inlet unsteadiness. An additional analysis was performed to show that the 
results obtained for the cascade of zero stagger angle was also applicable, 
with proper modification, to the cascades of arbitrary stagger angle. 


* The physical meaning of a complex number is attributed only to the real part, 
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Flow Around a Cascade of Nonstaggered Flat Plates 
for Inlet Flows with Translatory Oscillation 

Flow Configuration - Flow through the cascade shown in Figure 3 is considered. 
This cascade consists of flat plates of length, i, arranged with spacing, s, and 
zero stagger angle. A vertical translatory oscillation, w = w Exp (ivt), is 
superimposed on the steady through-flow, Uqq, in the direction of the flat plate. 
The steady inlet flow with a non- zero angle of attack is not considered here, 
because it can be easily added to the unsteady flow field. Since the flow is 
assumed to be two-dimensional, incompressible, and inviscid, the unsteady 
flow field through the cascade can be obtained by simultaneously solving 
Euler's equations of motion and the equation of continuity under the given inlet 
and boundary conditions. The potential flow field that satisfies the given inlet 
and boundary conditions at any instant agrees with the unsteady solution, as 
far as the velocity field is concerned. 

To determine the circulation of each airfoil, the Kutta- Joukowski theorem 
is assumed for the entire frequency range. The accuracy of this assumption 
is not clear in the higher frequency range because of the time lag involved in 
the formation of the boixndary layer around the airfoil. The circulation around 
each airfoil also fluctuates periodically, and this causes shedding of free 
vortices into the wake system. The wake of each airfoil, therefore, repre- 
sents a vortex sheet with varying intensity. The flow field induced by this 
vortex system must be taken into account. To simplify the analysis, wakes 
are assumed to be straight and parallel to the direction of through-flow, Uqq. 

The circulation of each airfoil, F, is divided into three parts, F = Fq 
+ Fj + Fj*. The flow with circulation Fq aroxind the airfoil induces a velocity 
distribution that tends to cancel the given inlet translatory oscillation, w 
= w Exp (ivt), on the airfoil. This solution represents a quasi-static flow when 
the effect of the vortex sheets is not taken into account. On the other hand, 
vortex sheets induce a downwash on each airfoil. The flow with circulation Fj 
induces a velocity distribution for canceling the non-uniform part of this down- 
wash on the airfoil, while the flow with circulation Fj ' induces a velocity dis- 
tribution for canceling the remaining uniform part. These conditions are in- 
dicated schematically in Figure 4. 

In this analysis, conformal transformation is used to obtain the potential 
flow solutions. 

Conformal Transformation - A cascade consisting of nonstaggered flat plates 
in the z = x + iy plane (physical plane) is transformed conformally into a circle 
of unit radius in the z* = x' + iy' plane by using the following mapping function 
(Reference 4); 


z 




R + z' 
R - z' 


+ In 


z' + 1/R \ 
z' - 1 / r ) 


( 1 - 1 ) 
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where R is a function of the solidity, cr = f /s, of the cascade; 


cr 


s 



R + 1 
R ^ 1 


(1-2) 


As shown in Figxire 5, the origin of the z-plane is located at the center 
of an airfoil, and the x-axis is chosen parallel to the airfoil. Therefore, 
z = ± i /2 describes the trailing and leading edge of the airfoil. These points 
transform into z* 1 in the z*-plane. Obviously, the flow aroimd the airfoil 
in the physical plane corresponds to the flow outside or inside of the unit circle 
in the z*-plane, and z = ± oo in the z-plane is transformed into z* = ± R or 
± 1/R in the z* -plane. 

Determination of - The flow field, which has the constant induced flow on 
the airfoil just opposite to the inlet translatory oscillation, w, can be deter- 
mined by the proper distribution of vortices in the z*-plane as illustrated in 
Figure 6. The solution, which includes no circulation around the airfoil, can 
be obtained in the z' -plane by locating vortices = - sw at z' = 1/R and 

at z* = - 1/R. On the other hand, the solution that produces circulation Fq 
around each airfoil is obtained in the z’ -plane by locating Tq/2 at z ' == ± 1/R 
and-Fo/2 at z' R. The assumed Kutta- Joukowski theorem requires that the 
total velocity induced by these vortices approach zero at the trailing edge (z' 

= 1), and this leads to the result 


ro = ^2 W Exp (ivt) (1-3) 

Determination of F| * - The circulation distribution in the wakes induces 
a downwash on each airfoil. This circulation distribution (Figure 4) can be 

00 

regarded as the sum of (1) Y-\^ on x = f /2 to oo and - F^ = - J dx at 

f/2 

X = + 00 , and (2) % at X = + 00 , The former distribution (1) induces the non- 
uniform part of the downwash, and the latter (2) induces the uniform part on 
the airfoil. The uniform part is calculated from the equation 

00 

w = Ty^/Zs - l/2s ^ y-\ff ^ 
i/2 

and the necessary circulation for canceling this downwash is determined by 
the same process that leads to equation (1-3); 


. 

+ 1 


2R 

R2 + 1 




(1-4) 
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Determination of F; - The solution for the non-uniform part of the downwash 
is described below. To simplify the analysis, the shedded vortices are 
assumed to be located on the straight lines starting from each trailing edge 
into the x- direction. First, the effect of the vortices, F^, at x = x is con- 
sidered. The uniform downwash of these vortices is subtracted by locating 
counter -vortices, - F^, at x = + oo as shown in Figure 7. This flow field can 
be obtained in the z* -plane by locating Fy^ at z* = x* and 1/R and - % at z* 

= 1/x’ and R, where x* is the mapped location of F^ calculated by the mapping 
function (equation (1-1)). 

To realize the Kutta condition at the trailing edge, circulatory flow around 
the airfoil is added to the above by locating Fj /2 at z* =± 1/R and -Fj /2 at z* 

= ± R, where Fj satisfies the following relation; 


R + 1 R - X 

+ 1 x' - 1 


(1-5) 


For the case of distributed vortices in wakes, Fj can be obtained by- 
integrating equation (1-5) along the wake; 


A 


= 2 


R + 1 
R2 + 1 


00 



i/2 


(1-6) 


In the above derivation, the circulation distribution of each vortex sheet 
is assumed to be identical. 

Unsteady Solution for Translatory Oscillation - According to Helmholtz's 
theorem regarding the persistence of vorticity, the total circulation of the 
whole system is assumed to be invariably equal to zero (Reference 5 and 6); 


r +^jT Y^dx = To + Fi + Ti* + = 0 (1-7) 

Putting equation (1-6) into eqioation (1-7), the following relation is obtained; 

f Ml frrr' 0 % 

f/2 

Since the inlet unsteady component is periodic with a frequency of f = v/2tt, 
the entire flow field also oscillates with the same frequency, f. It is also 
assumed that the vortices in the wakes move downstream with the velocity Uqq 
in the x-direction. Therefore, the circulation distribution, can be 

expressed in the form 
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( 1 - 9 ) 


= ~ g Exp [iv (t - VU^) ] 

where g is a constant complex number that describes the amplitude and the 
phase of the circulation oscillation in the wakes. Introducing the definition 


Fq + * = Gq Exp (ivt) 


( 1 - 10 ) 


wher^G^ is a constant complex number, the total circulation around the air- 
foil, P, can be written 


r = Exp (ivt) 


00 




r2 + 1 


R - x' 
x' - 1 


Exp (-iv x/Uqo) d 



( 1 - 11 ) 


Since the wake vorticity is produced by_the change of circulation around 
the airfoil, the increment of circulation, (dP/dt) dt, must be equal and opposite 
to the circulation in the wake between x = f /2 and x = i /2 + Uqo dt. Conse- 
quently, the relation dP/dt = - y^(f /2) Uqo is obtained; that is. 


dP . ^ . 

— = IV Exp (ivt) 


00 


C ^ R + 1 R - x' ^ . ,-rT s 1 /2 ttx\ 

^ s j 2 i2-n rr) 


f/2 


= (i /2) Uqo = - ~ g Exp (ivt) Exp (-iv f / 2 Uoo)Uoo 


The above treatment gives the correlation between G^ and g; 


00 


i/2 


( 1 - 12 ) 


. 2 irUoo 
- 1 


tUoo „ / . i \ 

Exp ( -IV -577- ) 

vs \ 2Uoo/ 


Using the relations derived from the mapping function (equation ( 1 - 1 )), 


7 



R - x' 

R + 1 

1 Exp 

er) 

- Exp (-ir<r) 

x' - 1 “ 

2 

^ Exp 

(^) 

- Exp (trcr) 


(R + 1)^ 2 

+ 1 1 + Exp (-■ir(r) 


equation (1-12) is written 


= Exp (-iTTCoa) 


[t^ 


C (cO, Q-) 


Exp (-ircr) 


t] 


(1-14) 


where 


CO = 


vs 


fs 


2 ^ 1 ^ ^ * reduced non-dimensional frequency. 


(1-15) 


C (o), cr) = Exp (iiTOJcr) 


00 

r 

1 

Exp “ Exp (-'[TO’) 

j 

e/z 

\l 

Exp - Exp (tto-) 






1 - t Exp (-Zttq-) 
1 ~ t 


- 1 


ioi - 1 , 
t dt 


(1-16) 


" L 


r(i + ito) 
r(i/2 + iw) 


n/tT 


r(i/2 + iw) 



ico. 


1/2 + ica; Exp (-2TTcr) 



The function, C (co, cr), represents the coupling effect between the wake and the 
airfoil and is expressed by using Gamma function, F (Reference 7 and 8), and 
Gausses hyper geometric function, F (Reference 9), as in equation (1-16). This 
coupling function is illustrated in Figure 8. In the case in which the solidity of 
the cascade is comparatively large, say o- g 1.0, the second term of C (co, cr) 
becomes very small and the function can be approximated only by the first 
term; 


C (cj, oo) = - 


A 

CO 



F(1 + ico) 
r(l/2 + ico) 


1 


(1-17) 


where C (0, oo) = In 4 = 1. 386. 

The integrated flow field can be obtained from the above result as follows: 

Upstream of Cascade - Since the total vortices included in the airfoils and 
wakes are invariably equal to zero (equation (1-7)), there is no induced velo- 
city far upstream of the cascade. Hence, the vertical velocity, Wj , far 


8 



upstream of the cascade is identical to the given inlet translator y oscillation; 


wj = w = w Exp (ivt) 


(1-18) 


Downstream of Cascade - The induced velocity in the y-direction at x = x 
downstream of the cascade, V 2 (x, y, t), is a ftmction of x, y and t. To 
simplify the analysis, however, this induced velocity is averaged over one 
spacing, s, thus giving 

s/2 

V 2 (x,t) = 1/s r V 2 (x, y, t) dy (1-19) 

-s/2 


As shown in Figure 9> an integration path enclosing the airfoil is set up 
extending from far upstream to x = x downstream. Using the definition of 
circulation, the following relation is obtained; 


- s V 2 = To + Fi + Fi ' + { dx 

i/2 


= Go Exp (ivt) + Fj + — g Exp (ivt) 

CO 


Exp (-iv tA) - Exp(-ir<r) 

' Uoo / 


( 1 - 20 ) 


From equation (1-4) and (1-7) the following equation is obtained; 


2R 




- (Fo + Fi + Fj') = - 


2R 


(R + 1) 


r (Ro + El) 


( 1 - 21 ) 


leading to 


— — 

Go Exp (ivt) = Fo - (Eo + Ej ) 


R^ + 1 ^ 
(R + 1)^ ° 


2R 

(R + 1)2 


El 


( 1 - 22 ) 


The first term of the above equation is given in equation (1-3), and the second 
term can be written as follows; 


A 


00 


= g Exp (ivt) 


I 

i/2 


R + 1 R - X' 
R2 + 1 x' - 1 



■ 2 C (co, (t) 

1 + Exp (-TT0') 

2C (cj.o-) _ 1 

_ 1 + Exp (-tto-) CO _ 


Go Exp (ivt) 


(1-23) 
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Substituting equations (1-3) and (1-23) into equation (1-22), C can be deter- 
mined: 


1 + io) 


Go = s [1 - Exp (-ircr)] 


2C (co,cr) 

1 + Exp (-TTcr) _ 
— w 


1 + io) C (o). O') 


(1-24) 


Using the relations obtained in equations ( 1-14), (1-23) and (1-24), each term 
of the right hand of equation (1-20) is related to the given inlet translatory 
oscillation, and thus the induced velocity, V 2 , becomes 


V 2 = - / ~ Exp [-iv(x - i/2)/Uoo] w Exp (ivt) (1-25) 

1 + 10) C (co, O’) 

The vertical velocity downstream of the cascade is calculated by adding the 
induced flow field to the given inlet translatory oscillation; 


W 2 = w Exp {ivt) + V 2 

= w Exp (ivt) |l - Exp [-iv (x - i/2)/Uoo]| (1-26) 

In the case in which the frequency of the inlet translatory oscillation approaches 
zero, 0 ) C{o),or) also approaches zero, and 


W 2 = W Exp (-TTO') (1-27) 

This result, of course, agrees with the solution of the steady flow. When the 
solidity of the cascade is not too small (or ^ 1.0), equation (1-26) can be 
approximated as follows; 

Wj = w Exp (ivt) I 1 - Exp [ -iv (x - !./2)/\3^]^ (1-28) 

This result is only good for the unsteady inlet flow with translatory oscillation. 

The solution obtained in the previous derivation represents a potential flow 
that is irrotational. In this case, the vortices are concentrated on wake lines 
starting from each trailing edge, and move downstream with velocity Uoo , as 
shown in Figure 10. When the cascade follows another cascade, the flow down- 
stream of the first cascade becomes the flow upstream of the second. The 
real outlet flow from a cascade, therefore, imposes a very complicated inlet 
condition to the following cascade. 
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To avoid these difficulties, the real outlet flow is replaced by a simplified 
outlet flow. The simplification is made by redistributing vortices of dx 
on wakes (length = dx) uniformly over the region, s by dx, as shown in Figure 
10. This simplified flow is no longer irrotational; however, it has a vertical 
velocity that is independent of ordinate y. The mean vertical induced velocity, 
V 2 , of equation (1-25) also describes the vertical velocity of the simplified 
flow, which is nothing but the sinusoidal gust oscillation. 

Flow Around a Cascade of Nonstaggered Flat Plates 
for Inlet Flows with Sinusoidal Gust Oscillation 

Flow Configuration - The inlet flow has a vertical sinusoidal gust oscillation, 
w = w Exp [iv(t - x/Uoo)], propagating with the main through- flow Uoo. J’he^ 
unsteady circulation, f , around the airfoil is divided into thre^parts, F = Fq 
+ Fj + Fj', as shown in Figure 11. The flow with circulation Fq induces the 
down wash which cancels the given inlet sinusoidal gust on the airfoil so that 
the surface of the airfoil becomes a part of the instantaneous stream line. 

The flow with circulation Fj and Fj * cancels the induced flow by vortices in the 
wakes and, therefore, has the same relations as in the case of the translatory 
oscillation. 

Determination of Fq - As shown in Figure 12, the flow in the physical plane is 
conformally transformed into a unit circle in the z' -plane using the mapping 
function (equation (1-1)), To obtain a sinusoidal downwash opposite to the gust 
on the airfoil, a more complicated distribution of vortices in the z' -plane is 
necessary than in the case of a uniform downwash, A circulation distribution, 
V(0), in the form of 

00 

y(0) = ^ Aj^ cos n0 (1-29) 

n=l 


is located on the unit circle. The Fourier coefficients, Aj , A 2 , . . . , will be 
determined so that the induced velocity on the airfoil is opposite to the given 

Ztt 

inlet gust. Since the sum of the distributed circulation, ^ 7(0) d0, is equal to 

0 

zero, the circulation distribution produces no circulatory flow around the air- 
foil as a whole. 

An additional circulation is superimposed on the flow by locating Fq/ 2 at 
z* =± 1 /R and - Fq/ 2 at z* = ± R to satisfy the Kutta- Joukowski theorem at the 
trailing edge. 

The velocity (u*, v*) induced by the circiolation distribution Y(0) at z* = 1 is 


(u* 


iv») 


z*=l 


00 

i V r - ^n cos n0 
2 tt Z J 1 - Exp (i0) 
n=l 0 
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Since the integrand becomes infinity at z* = 1 (0 = 0), this integration is 
carried out along the path shown in Figure 13 and gives 


(u* 



00 



n = l 


lim 
6 ^ 0 


00 


27T-€ 


i Z I ^ 


cos n0 


Exp (i0) 


d0 


1 = 1 6 


the above = 0 


A r 

n J 1 - Exp (i 


0 ) 


00 


2 ^ 


. n=l 
= 1 T" 


(1-30) 


Hence, the induced velocity at z* = 1, 


00 

\ 


u* = 0, V* = - 


Z 

n = l 


must be canceled by the proper circulatory flow, and this necessarily leads 
to the result 




00 00 
^ A = - 7T sech (ttct/Z) ^ 


R'" + 1 n 

n = l 


n=l 


(1-31) 


The velocity induced by the circulation 7(0) at an arbitrary point on the unit 
circle, z’ = Exp (i0Q ), is written as follows; 


(u* - iv')_^, 


lim 

z' € -► 0 


00 00 +277-e 


i I 


An cos n0 


Exp (i0o ) - Exp (i0) 


d0 


00 


n = l 00 + ^ 
0o+€ 


t; Z ^n=os„e. J 


d0 


Exp (i0o) - Exp (i0) 


n = l 


0o“6 


00 


= -| Exp (-i0o) ^ Exp (-inOg) 


(1-32) 


n = l 


The velocity field (u, v) in the z-plane can be calculated from the velocity in 
the z* -plane using the relation 


(u - iv) = (u* - iv' ) 


dz^ 

dz 


(1-33) 
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dz ^ 

where — is derived from the mapping function (equation (1-1)); 
dz 


/dz'\ _ ^ - Exp(2ien)1 rExp(2ien) - 1/RM 

^<*“''z'=Exp (ie„) ° ® ^ 


Substituting equations (1-32) and (1-34) into equation (1-33), the expression 
is obtained for the induced velocity in the z-plane; 


u - iv 


00 

^ (R + 1/R)^ - 4 cos^ 0 Y 
s 8 (R + 1/R) sin 0 L 

n=l 


Aj^ Exp (-in0) 


(1-35) 


or 


V 


2ir 

s 


(R + 1/R)^ - 4cos^0 
8 (R + 1/R) sin0 




n=l 


sin n6 


(1-36) 


at a point on the airfoil corresponding to z* = Exp (i0) in the z* -plane. In the 
previous derivation, all of the Fourier coefficients are confined to real num- 
bers, If not, the horizontal induced velocity, u' , at z* =1 can take a certain 
finite value that cannot be canceled by adding circulatory flow around the air- 
foil, and the Kutta- Joifkowski condition cannot be satisfied. 

Next, the complex Fourier coefficients, A^, are introduced and replace 
the real coefficients, A^, in equation (1-36) with the real part of A^^ Exp (ivt); 


V 


^ (R + 1/R)^ - 4 cos^e 
s 8 (R + 1/R) sin 0 




Aj^ Exp (ivt) sin n0 


(1-37) 


Since the circulatory flow around the airfoil with circulation Fq induces 
no vertical velocity component on the airfoil, equation (1-37) describes the 
total vertical component. To satisfy the boundary conditions on the airfoil, 
the vertical induced velocity must be equal and opposite to the given inlet 
gust. Thus, 


V + w = V + w Exp (ivt) Exp (-iv x/Uqo? = 0 


(1-38) 


Substituting equation (1-37) into equation (1-38), the following equation is 
obtained; 


I 


Aj^ sin n0 


s 


8 (R + 1/R) sin 0 
(R + 1/R)^ - 4 cos^ 0 


w 




(1-39) 
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Since the above function, f(0), is symmetrical about 0 = tt, the function can 
be expressed only by a sine series. The Fourier coefficients, Aj , A 2 , . . . . , 
are determined by 


2tt 


An = 1 /tt ^ f(0) sin n0 d0 


0 

Hence, the sum of the Fourier coefficients is 
00 00 

^ A„ = 1/, ^ r £(0) sin n0 d0 


n=l 


n=l 0 


(1-40) 


(1-41) 


Using integration by parts, the following equation is obtained; 

2tt 

= ./,]• £.,e,Ji£^^d0 


n=l 


n=l 


the above = 0 


n=l 


The relation, 
00 


Zj n 


n0 


= - In (2 sin 0/2), for 0 < 0 < 2 tt 


n = l 


leads further to the following expression; 

00 2 tt 

^ = - 1/tr J f (0) In (2 sin 0/2) d0 


n=l 


2tt 


= - l/iT I f(0) In (2 sin 0/2) | + l/2irj f(0) cot (0/2)d0 


the above = 0 


0 


This integral will be finally conducted in the z-plane after being transformed 
by the following relations; 


s 4 (R + 1/R) sin 0 

^0 = dx. 


2tt (R + l/R)^-4cos^0 


cot 0/2 = Exp (iT(r/2) 


Exp - Exp (-'iro') 

Exp (tto-) - Exp 
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The result is 




- “ Exp {-tt(t/2) [ Exp (iT(r) - Exp (-ucr)] 

TT 


(1-42) 

Exp (-iTTCo(r) S (go, (t) w 


where 


_ Exp (iTcr) 

'O’ Exp (ttct) - Exp (-TTcr) 


Exp (i-rrcocr) 


i/2 

x=-i/2 


Exp - Exp (-ircr) 

Exp (iror) - Exp 


1 + icj 


dt, 



B = l/[Exp (2tro-) - 1] 


The sinusoidal gust function, S(cj,cr), can be calculated by expanding the 
integrand into a power series of t, and is written as follows; 


S(co,o-)=^F[l/2, 1+icj, 2; 

_ tanh (ttcjo) E( 1 + ico) 

” oj IXl/2 + io)) 


1 - Exp (-2ttct)] 
F[l/2, 1 + ico. 


1/2 + ico; Exp (-2ircr)] 


(1-43) 


Exp (■TTcr) Exp (2iTrcacr) 
n/tt” (1 - 2iu)) 


r(i/2 + ico) 
F(1 + icj) 


F[3/2, 


1 - ioo, 3/2 - ico; Exp (-2'rrcr)] 


This function is illustrated in Figure 14. As seen from equation (1-43), 

S(co, O’) oscillates slightly with the reduced frequency co, since it contains a 
periodic term, Exp (Zincoa). When the solidity of the cascade is not too small, 
this oscillation is not essential. In the case in which the solidity is compara- 
tively large, this function can be approximated by 


s (w.oo) F(l/2, 


1 + ica , 2 ; 1 ) 


tanh (ttco) F(1 + ico) 
^ CO rtl/2+i(o) 


= complex conjugate of ^ ^ ^ (co, «) 


(1-44) 
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In the limiting case of co = 0, the function becomes 


S (0,(t) = 


1 

1 + Exp (-TTcr) 


(1-45) 


By substituting the real part of 


I 


Exp (ivt) into equation (1-31), the 


n = l 


circiilation around the airfoil is found to be 


Pq = 2 s [1 - Exp (-TTcr)] Exp (-iTtcoij) S (co, cr) w Exp (ivt) (1-46) 

Unsteady Solution fo r Sinusoidal Gust Oscillation - As illustrated in Figure 11, 
the total circulation, F, aroiind the airfoil consists of Fq, Fj and Fj'. Since 
the total circulation of the system remains invariably zero. 


oo 


i/2 



0 


An analysis similar to that leading to equation (1-14) yields the expression 


= Exp (-iTTCocr) 


2C(co,o-) 

1 + Exp (-7r(7) 


CO 


where 

Gq Exp (ivt) = fo + El' , 

27T 

^ Exp [iv (t - VU^) ] 

The integrated flow field for a given inlet sinusoidal gust oscillation, 
w = w Exp (ivt) Exp (-iv x/Uqo), can be obtained as follows: 

Upstream of Cascade - Since the total vortices included in the airfoils and 
wakes are equal to zero, there is no induced velocity at a sufficient distance 
upstream of the cascade. Therefore, the vertical velocity, Wj , far upstream 
is equal to the given inlet sinusoidal gust oscillation; 


Wi = w = w Exp (ivt) Exp (-iv x/Uoo) 


(1-47) 
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Downstream of Cascade - The induced velocity, V£ , at x = x downstream of 
the cascade is determined in a manner similar to the case of translatory 
oscillation; 


- s V 2 = To + A + Fi' + J dx 

i/2 


(1-48) 


= Gq Exp (ivt) + Fj 



g Exp (ivt) [Exp (-iv x/Uoo) - Exp (-irtr)] 


Further, 


Go Exp (ivt) = 


+ 1 

(R + 1)=^ 


ro 


2R 

(R + 1)2 


A 


(1-49) 


where 


Fq = 2 s [1 - Exp (-TTO')] Exp (-iiTCocr) S (co, a) w Exp (ivt). 


ri 


“ 2 C ((U),cr) 

1 + Exp (-TTcr) 

2 C (co, (r) 

1 + Exp (-TTcr) 00 


Go 


Exp (ivt) 


From the above. 


1 + i^ , ^ — ( ) 

s [1 - Exp (- 2ttct)] Exp (-iTTCocr) S (co, cr) — - ~ w (1-50) 

■' 1 + ICO C (co, cr) 


2 C (co, cr) 


Go = 


Each term of equation (1-48) can be expressed in terms of the given inlet 
sinusoidal gust oscillation, and the induced velocity, V 2 , becomes 


V 2 = - [1 - Exp (-27T(r)] — ‘v w Exp (ivt) Exp (-iv x/Uqo) (1-51) 

i + 100 O ^00 , or j 

where V 2 represents the mean value on one spacing, s. Finally, the vertical 
velocity downstream of the cascade is 


W 2 = w Exp (ivt) Exp (-iv x/Uoo) + ^2 


(1-52) 


= w Exp (ivt) Exp(-iv x/Uoo) |l - [1 - Exp (-Zircr)] ^ + f <r)} 
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For a comparatively large value of solidity (<r g 1.0), equation (1-52) can be 
approximated by 


W 2 = w Exp (ivt) Exp (-iv x/Uoo) 1 


tanh (ttco) 

TTO) 


(1-53) 


When CO = 0, coC (0, cr) = 0 and S (0, cr) = [ 1 + Exp (-ucr)]*^, thus W 2 = w Exp (-ttct). 
This result agrees with the solution for steady flow. The above result gives 
the unsteady flow solution only for an inlet flow with sinusoidal gust oscillation. 
If the inlet flow has combined translatory and sinusoidal gust oscillation, the 
solution can be obtained by superimposing both flow fields. 

When the inlet flow oscillates periodically, the lift acting on the airfoil 
also oscillates periodically. This unsteady lift is calculated for both trans- 
latory and sinusoidal gust oscillation, and is given in the Appendix of 
Reference 10, 

Solution for a Cascade with an Arbitrary Stagger Angle 

Flow Approximation aiui J.ts ^c cujr^cy - The dynamic property of a cascade 
having an arbitrary stagger angle are analyzed. Since the exact solution for 
this case can be expressed only with impractical mathematical complexity, 
an approximate method is adopted. Flow through a staggered cascade of 
flat plates, as shown in Figure 15 (a), is considered. The inlet flow is a 
superposition of through-flow, Uqo, in the direction of the flat plate (x-axis) 
and a translatory or sinusoidal gust oscillation in the direction of the cascade. 
First, the flow fields through nonstaggered and staggered cascades having the 
same solidity are compared (Figure 16). Both inlet flows are assiomed to be 
steady and have the same through-flow, Uqoj nnd the same w in the cascade 
direction. 

The influence coefficient, k, of the staggered cascade (Reference 4) is 

Fo- 4s R cos Qt /I 

“ To- = 0 “ iri +1 cos X 


where R and 0,^ are characteristic numbers included in the mapping func- 
tion that transforms a staggered cascade into a unit circle in the z’ -plane, as 
illustrated in Figure 15 (b). R and 0rp are determined by solving the 
following equations simultaneously; 


cr 


i/s = 1 



cos \ In 


r2 


+ 2R cos 0T 
- 2R cos 0x 

+ 2 sin X tan“ 


+ 1 
+ 1 

1 sin 0T 

- 1 


■) 


tan 0T = 


R^ - 1 
R2 + 1 


tan X 


(1-55) 
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The ratio of the circulation of staggered and nonstagger ed cascades is written 


_ ^ ^/£ . or .- = . 0 ^g- .. = . Q .. ^ _!< — ^ . 03 . 0 .. ■ (1.56) 

^\=0 = o/^\,ar = 0 ^\,o- =0 = 0 ^\,cr=0 

Since ^ - 0 /^\ (r = 0 indicates the circulation ratio of the isolated airfoil, 
this ratio is proportional to the vertical velocity component in the y-direction, 
and thus. 


^ cr = 0 
cr = 0 


w cos X 
w 


= cos \ 


(1-57) 


From equation (1-54), (1-56) and (1-57), 

R 

F p2 1 1 COS 0'p 

-= = — — = function of (<r, \) (1-58) 

= o 

Figure 17 shows the numerical values of equation (1-58), Note that the 
circulation ratio remains close to unity if the solidity, (r , does not become 
too small. From the definition, Fo/(Fo)^_q and Fi*/(Fi')^_q are exactly the 
same as F/F^_q of equation (1-58) and are approximately unity. Through a 
similar comparison it can also be concluded that Fo/(Fo)\-q and Fj/(F|))^-q 
remain close to unity if the solidity, cr , does not become too small. The 
above comparison shows that the result obtained for a nonstaggered cascade 
can also be applied to a staggered cascade with reasonable accuracy, provided 
that the unsteady flow component in the direction of cascade is the same for 
both cascades. 

This approximation is, of course, not sufficiently accurate for determin- 
ing the steady flow solution. However, since we are mainly interested in the 
dynamic characteristics of cascades, this approximate method seems to be 
reasonable. 

Results 

Translatory Oscillation Case 

At a sufficient distance upstream of the cascade; 

Uoo in the x-direction, 

Wj = w Exp (ivt) in the cascade direction. 

Downstream of the cascade; 

Uqq in the x-direction. 
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(1-59) 


W2 


= w Exp (ivt) s 1 


ri - Exp (-rro-)l _ f . v 
i + icoc(co,^) 


(x-y tan \ ■ 


f>]} 


in the cascade direction, where 


vs 

ZirUoo 


vs 


COS \ 


Sinusoidal Gust Oscillation Case 

At a sufficient distance upstream of the cascade; 
Uoo in the x- direction. 


wi = w Exp (ivt) Exp 


Uo 


(x-y tan X) 


in the cascade direction. 


Downstream of the Cascade 
Uoo in the x-direction. 


W2 = w Exp (ivt) Exp 



(x-y tan 



[1 - Exp (-Zttct)] S (co, (t) 
1 + icu C (co, 0 *) 


} 


in the cascade direction 


(1-60) 
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SECTION II. ANALYTICAL STUDY 
OF DYNAMIC CHARACTERISTICS OF TURBOPUMPS 


Using the results obtained in the previous Section, a theoretical analysis 
was performed to determine the relationship between pressure rise and flow 
rate through a turbopump when the flow rate fluctuates periodically aroimdits 
mean value at a constant frequency and while the rotational speed of the turbo- 
pump is constant. 

There is a wide variety of principles of turbopump construction ranging 
from the single-stage axial flow pump to the multi-stage centrifugal piamp. 

Since the working principle of the axial flow pump is based purely on the per- 
formance of a linear cascade, the axial flow pump offers the most fundamen- 
tal case for theoretical application. Therefore, a method for determining the 
dynamic characteristics of a single-stage axial flow pump is described first. 
This method serves as an example of the theoretical approach for determining 
the dynamic characteristics of an arbitrary turbopump. A simple method is 
also derived for the case when approximate dynamic characteristics are 
sufficient. This method determines the dynamic characteristics without enter- 
ing into the details of pump design. 

In these analyses, the flow was assumed to be incompressible and inviscid. 
The effect of eventual cavitation in the pumps is not considered. 

Dynamic Characteristics of a Single-Stage Axial Flow Pump 

Pressure Difference Due to Oscillating Through- F low in Conduit - The analy- 
sis was conducted on the axial flow pump illustrated in Figure 18. This pump 
has only one stage consisting of a rotating blade row (rotor) and a following 
stationary blade row (stator). The rotor is driven by a prime mover at a con- 
stant rotational speed, n. The flow rate through the pump, Q, is assumed to 
fluctuate periodically with a frequency, f, aroimd its mean value, Qq. Assioming 
also the sinusoidal fluctuation, the instantaneous flow rate is 

Q(t) = Qq + AQ Exp (ivt), (2-1) 


V = Zivf 


where the physical meaning of a complex number is attributed only to its real 
part. 

In the actual pump, the flow pattern around the rotor and stator varies 
with radius, r. For the calculation of dynamic characteristics, however, a 
representative radius, is chosen for the analysis as shown below; 


r 


m 



+ rh^ 
2 


( 2 - 2 ) 
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where and are outer and inner radii of annular passage. This repre- 
sentative radius has been used often for the analysis of steady- state charac- 
teristics and has proved to give good results (Reference 11). 

By developing the cylindrical flow surface of radius r^^ into a plane, a 
two-dimensional flow is obtained as shown in Figure 19, where the rotor cas- 
cade moves straight in the cascade direction with a constant velocity, u = 
Zirnr^^. 

Since an incompressible, two-dimensional flow is assumed, the axial 
velocity, remains constant upstream and downstream of the rotor and 

stator. 

The flow pattern of an axial flow pump is in most cases the free vortex 
type, especially when the pump is not provided with inlet guide vanes as 
illustrated in Figure 18. 

The free vortex type flow pattern assures the constant axial velocity over 
the entire radius, and thus 


Va(‘) = , 2, = Exp (ivt) 


wher e 

V 


Q 


Q- 


AQ 




(2-3) 


The mean velocity triangle at each station is illustrated in Figure 19. 

The pressure difference between the delivery and suction ports of the 
pump, (p^ - Pg), oscillates periodically corresponding to the fluctuating flow 
rate. As the first step, the pressure difference will be calculated when there 
are no cascades between suction and delivery ports. This case represents an 
oscillating flow in a conduit with varying sectional area as shown in Figure 20. 
If only one-dimensional flow in the conduit is assumed, the velocity in the con- 
duit, V(X, t) is 


V(X, t) 


Q(t) Qo 

A(X) A(X) 


+ 


AQ 

A(X) 


Exp (ivt) 


(2-4) 


where A(X) is the sectional area at distance X downstream of the suction 
reference point. 

The pressure gradient, 9p/9X, in the conduit can be calculated by substi- 
tuting equation (2-4) into the equation of motion for a non- viscous, incom- 
pressible fluid 

I iR - ^ V — 

" P 9X “ 8t ax 
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(2-5) 



In the usual design practice, the change rate of sectional area, 9A(X)/9X, 
is made as small as possible to reduce fluid losses, and this results in a 
small change rate of through- flow velocity, 9V/9X. 

From the above, the first term of equation (2-5) is predominantly larger 
than the second term, and consequently 


9p 9V 

8x " ■ ^ at ~ 


. A Q 


Exp (ivt) 


By integrating the above pressure 
delivery port. 


gradient from the suction port to the 


L 

Pd - Ps = - dX 

X =0 


= - piv AV^ Lgq Exp (ivt) 


(2-6) 


where 

Aq = •ir(r^^ - rj^^) : area of annular passage, 

L 

Eeq = I dX : equivalent length. 

X = 0 

The above pressure difference occurs independently of pumping action 
and is regarded as an inherent property of the conduit. 

The dimensionless pressure difference, h^, indicating the conduit effect 
is thus 


he = 




2iv ^ Exp (ivt) 

u u 

L 0 Q AQ 

2ico ^ -pr — 4> Exp (ivt) 

^m ^ o 


(2-7) 


where 


CO 


vr 


XQ- 


U 


— : reduced frequency of pump, 
n 


C|) = 


flow coefficient. 


u 
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As seen from equation (2-7), always varies with 90® delay to the 
fluctuation of flow rate throughout the entire frequency range, and the magni- 
tude is proportional to the frequency, f, at a constant fluctuation ratio, AQ/Q^. 

Change of Pressure Rise Through Rotor - The rotor cascade consists of 
arbitrary airfoils arranged with arbitrary spacing and stagger angle as 
indicated by a broken line in Figure 21. This cascade is then substituted by 
a cascade of flat plates that have the same chord length as that of the airfoil 
and are located parallel to the zero-lift inlet direction of the original airfoils 
in the cascade. 

The absolute inlet velocity into the rotor has no component in the cascade 
direction and has only the axial component = ^ao ^^a which 

oscillates periodically. 

The relative inlet velocity, W, into the rotor is obtained by superimpos- 
ing axial velocity on the constant peripheral velocity, u, vectorially. 

It is now assumed that the flow field corresponding to the mean steady 
inlet flow, W 2 o> has been solved already by an appropriate method (Reference 
12 and 13), and that the mean pressure rise through the rotor has been deter- 
mined from this flow field. The deviation of pressure rise from its mean 
value is determined below. 

The unsteady axial component, Exp (ivt), upstream of rotor is divided 

into two components, AV^ sec Exp (ivt), in the direction parallel to the flat 
plate, and wj = - AV^^^ tan Xp^ Exp (ivt), in the cascade direction where Xp is 
the stagger angle of the simplified rotor cascade. 

The component, AV^^ sec Exp (ivt), which is parallel to the zero-lift 
direction of the original cascade, induces no circulation around airfoils and 
thus passes through without causing any change of flow field. On the other 
hand, Wj in the cascade direction represents a translatory oscillation of inlet 
flow and induces a fluctuating circulation around the airfoils. 

The solution for this case is already obtained in Section I as follows: For 

upstream of rotor, 

W| = w Exp (ivt) 
wher e 

w = - AV^ tan Xp (2-8) 

For downstream of rotor, 

W2 = Wj + V2 (2-9) 
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where 


V 2 tan \R iVicfR^C(™lR) Exp(ivt)Exp [-iv (X - ^ cos \r)/V^^J . 

_ _ cos \r 

~ ZirUoo " Nr ()> “ 

In the above relation, is the spacing of the rotor cascade and Nj^ is the 
number of rotor blades. The coordinate systems, (x, y) and (X, Y), are 
illustrated in Figure 21. 

The flow field downstream of the rotor cascade is obtained by superimpos- 
ing sec Exp (ivt) in the x-direction and W 2 of equation (2-9) in the Y- 

direction. The superposition of sec Xj^ Exp (ivt) and the first term of 

equation (2-9) produces a simple oscillating flow, AV^ Exp (ivt), in the axial 
direction that is equal to the given inlet unsteady condition. The disturbance 
induced by the rotor cascade exists only in the second term of equation 
(2-9). The flow field downstream of the rotor is, therefore, the superposition 
of the mean flow, W 20 (steady); axial oscillation, AV^^ Exp (ivt); and the induced 
disturbance, V 2 , in the cascade direction as illustrated in Figure 21. The 
associated pressure change is calcxilated from this flow field. 

The basic equation for determining the pressure of incompressible, potential 
flow is 

£ + i5 + S2+^v2= F(t) (2-10) 

p at 2 

wher e 

$ velocity potential 

^2 potential due to gravity 

V resultant velocity of fluid 

F arbitrary function to fit initial condition 

The gravitational potential, Q, produces merely the hydrostatic pressure field. 
Therefore, the potential, is omitted hereafter by considering the difference 
of real pressure and hydrostatic pressure. 

To calculate the pressure by equation (2-10), it is necessary to determine 
the so-called "impulsive pressure, p(a$/at)." Since the velocity potential of 
the mean steady flow is also steady with respect to time and does not contribute 
to impulsive pressure, only the velocity potential of unsteady flow is significant. 

The unsteady velocity potential is divided into three parts; is due to 
through flow, $2 due to vertical translatory flow^and #3 is due to induced flow 
as shown in Figure 22 (a), (b) and (c). 
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The super imposition of through- flow and vertical translatory flow gives 
the oscillating flow in the axial direction. Using the X, Y coordinate system 
shown in Figure 21, 


u = (ivt) in X-direction 

V = 0 in Y-direction 
and thus, 

+ $2 = ^ Exp (ivt) 

Differentiating the velocity potential as to time, the impulsive pressure is 


a($7 + $.) 
8 t 


piv AV 3 ^ X Exp (ivt) 


This impulsive pressure results in the same pressure difference given 
by equation ( 2 - 6 ), and was already taken into account as a part of h^ in equa- 
tion (2-7). Therefore, the contribution of and $2 ^o the impulsive pressure 
need not be considered again. 

The flow corresponding to $3 is induced by vortices in the airfoils as well 
as in the wakes and has the velocity and pressure distribution that are a func- 
tion of X, Y, and t downstream of the cascade. However, as far as this real 
downstream flow is approximated by the simplified flow as shown in Figure 
22 (d), it is obvious from the equation of motion that the pressure, P 2 , is 
instantaneously constant anywhere downstream of the cascade. Thus, the 
pressure downstream of the cascade is equal to the mean pressure at the exit 
line of the cascade at any instant. 

The detailed analysis on the effect of $3 on the pressure at the exit line 
of the cascade (Reference 14) leads to the conclusion that the impulsive pres- 
sure, p( 8$3 / 8 t), produces a local pres sure gradient between two adjacent blades, 
but its mean value over one spacing becomes zero. Since only the mean 
pressure over a spacing is of interest, the contribution of ^3 to the down- 
stream pressure, can be omitted from consideration. 

It has become obvious that the impulsive pressure term, as well as the 
gravitational potential term, can be dropped from equation ( 2 - 10 ). 

The pressure is finally calculated by 


£ 

P 


^ - F(t) 


( 2 - 11 ) 


From the velocity triangle in Figure 21, the instantaneous velocity head 
upstream of the rotor is 
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1/2 {[V^o + AV^ Exp (ivt)]2 + 


where the subscript u denotes the component in the cascade direction. 

On the other hand, the velocity head at the exit of the cascade is from 
equation (2-9) 


>« j[vao 


+ AVg^ Exp (ivt) 


2 


+ 




1 - Exp (-ttctr) 

1 + icoj^ C(coj^, CTj^) 


Exp (ivt) 



From equation (2-11), the deviation of pressure rise from its time-mean 
value, App^, is 


^PR 

P 


= 1/2 


K . o ' 


w. 


u 


20 


+ AVg^ tan Xp^ 


ri - Exp (-iro-R)l 

1 + C(wj^, o-j^) 


Exp (ivt) 



- 1/2 (W. 


aio 


- W.. 


^2 0 


W, . AV^ ^ Exp (ivt) (2-12) 

^ ® 1 + i^p^ C ^ 

In the above, AVg^ is assumed to be small and the term of AVg^^ is neglected. 
The dimensionless pressure change is then 


hp = 


^PR 


- 2 tan Xp 


1 - Exp (-iro-p) Q ^ 

1 + ^^R ^ ^^R • ^R^ ^ ^o 


<j> Exp (ivt) 


(2-13) 


Change of Pr essure Rise Through Stator - The inlet flow to the stator contains 
two kinds of inlet unsteadiness. The oscillating axial component, AV^^^ Exp (ivt), 
produces a translator y inlet oscillation, Wj = - AVg^ tan Xg Exp (ivt), to the 
stator similar to the case of the rotor. Besides this, the distiirbance velocity 
induced by the rotor generates a simulated sinusoidal gust oscillation upstream 
of the stator and makes the flow through the stator more complicated than the 
flow through the rotor. 

The flow downstream of the stator includes two kinds of distrubance, one 
resulting from translatory inlet oscillation, and the other from sinusoidal gust 
oscillation. 

The flow caused by the oscillation of axial flow is solved in exactly the 
same way as in the case of the rotor. The flow downstream of the stator is 
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represented by 


AVa Exp (ivt) in the X-direction, 

V 2 = AV» tan Xg (ivt) Exp j^-iv(X - ^ cos Xg)/V^„] 


in the Y-direction 

where coordinate (X, Y) is illustrated in Figure 23. 

On the other hand, the flow upstream of the stator contains a sinusoidal 
gust expressed generally as 


Wj = w Exp (ivt) Exp (-ivX/Vg^Q) 

where w must be specifically so determined that the disturbance induced by 
the rotor can be properly transmitted into the coordinate system of the stator. 
From the condition that the above equation at X = - [d + f g/2(cos \g)] coincides 
with equation (2-9) at X =fj^/2(cos \^), the inlet sinusoidal gust to the stator 
has the form 


Wi = - AVg^ tan 


1 - Exp (-•(^(tr) 

1 +iwR C(coj^, o-p^) 


Exp 


Is 


iv(d + cos Xg)/V 


ao 


(2-14) 


Exp (ivt) Exp (-ivX/V^^) 


The solution for a given inlet sinusoidal gust is solved in the previous 
Section in which the flow, w^ , downstream of the cascade is known 


w. 


+7^ =5. |l - [1 - Exp (.Zxcg)] 1 fts C(os, vs)} 


where 

vss cos Xc; 
S ZttUoo Ng<j) 


and 


sg is the spacing of the stator cascade, 

Ng is the number of the stator blades. 

The first term of equation (2-15) is a sinusoidal gust oscillation, which 
merely passes through without any change, while the second term indicates 
the induced disturbance because of the inlet gust flow. The entire flow field is 
obtained by superimposing the steady mean flow, the translatory oscillation 
flow, and the sinusoidal gust oscillation flow as illustrated in Figure 23. 
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For the calcxilation of pressTire, it is necessary to consider the effect of 
impulsive pressure, p( 8 $/ 8 t). So far as the translatory oscillating inlet flow 
into the stator is concerned, the effect of impulsive pressure is restricted to 
the generation of pressure gradients within a spacing as in the case of the rotor, 
and thus can be omitted from the consideration. 

Besides the translatory oscillation, the flow around the stator is accom- 
panied by the sinusoidal gust oscillation, which is no longer irrotational flow 
and for which the velocity potential does not exist. This circumstance makes 
the determination of pressure very difficult, since the fundamental relation 
for pressure, equation (2-10), does not hold for rotational flow. To avoid this 
diffic\ilty, the real sinusoidal gust inlet flow is replaced by the modified flow 
field as illustrated in Figure 24 (a). 

In the modified flow, the inlet flow has no gust oscillation and is conse- 
quently irrotational. Instead of this steady inlet flow, the airfoil itself 
oscillates with the velocity - Wj in the cascade direction, so that the airfoil 
encoimters the same relative velocity to the fluid. Even through this exchange 
of inlet and bovindary conditions, the same circulation distribution will be 
obtained in the airfoil and in the wake as that of the original flow, provided 
|wi I « This replacement, however, has an essential advantage because 

there exists a velocity potential. Therefore, the effect of impulsive pressure 
will be calculated from the modified flow field. 

The unsteady velocity potential of modified flow is divided into two parts, 

$4 and $ 5 . The potential $4 is generated by the unsteady vortices in the air- 
foil and wake that are conditioned by the oscillating motion of the airfoil. 

An analysis similar to that leading to the result p( 9 $ 3 / 8 t) has no effect on 
an average pressure of rotor cascade, gives the conclusion that the impulsive 
pressure, p( 9 $ 4 / 9 t), generates merely a pressure gradient between two 
adjacent blades, and has no effect on the average pressure over one spacing. 
Since we are not interested in the local pressure but rather in the average 
pressure, the effect of impulsive pressure, p( 9 § 4 / 9 t), can be dropped from 
equation ( 2 - 10 ). 

Besides the tinsteady circulation, each airfoil has a steady, distributed 
circulation, Vo(^)> over its airfoil. The total steady circulation Tq of the 
airfoil is the sum of distributed circulation, and thus, 

“ cos \ 

^o= y Yo(X)dX = (V^^^ - 

X = - ^ cos\ 

In the case when the airfoil does not move, this steady circulation causes 
no change of velocity potential as to time and produces no impulsive pressure. 
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In the modified flow, however, the airfoil is in an oscillating motion, 
and this causes an associated time- dependent velocity potential field, § 5 . 

Now the velocity potential can be calculated of a strip of cascade with 
width dX, which is moving with velocity -Wj in the cascade direction, cf. 
'Figure 24(b). The time -differential of the velocity potential at point 1 
sufficiently upstream of the strip is 

/BY\ Yo dX ^ 

(.if). = (if), (if). = ^ 

Since only the relative motion between point 1 and the strip is essential, 
the point 1 is supposed to move with velocity dY/dt = Wj in the above deriva- 
tion, instead of considering the strip as moving with -Wj . 

From the same analysis, the time -differential of the velocity potential 
at point 2 sxifficiently downstream of the strip is 



Yq dX 
2 s 


Wj 


From the above, the difference of time -differential across the cascade 
becomes 


The integration of equation (2-16) over the airfoil gives the total difference of 
time -differential due to total circulation on the airfoil; thus. 


— cos \ 

(If^). -(^). =F I v„(X)«. (x,.)dx 

X = - COS \ 

“ cos \ 

= -J Exp (ivt) J Yo (X) Exp (-iv X/V^q) dX (2-17) 

- “ COS \ 


To carry on the above integration, the distribution of Yq (^) must be given, 
which is different in each case according to the geometrical shape of the airfoil. 
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It is assumed here, however, that the circulation distribution changes linearly 
from its maximum at the leading edge to zero at the trailing edge. Thus, 


7o(X) = ^ (l - . (2-18 

^ 1 cos \ \ i COS \/ 

Substituting equation (2-18) into equation (2-17) and integrating over the 
airfoil, we obtain 

\ dtj^ 


vs cos \ 

" 2ttUoo “ 2tt 

This result should be considered for the calculation of pressure around 
the stator. 

The pressiore downstream, of the exit edge of the stator is consider ed to be 
constant for the same reason as for the rotor, and the pressure difference 
upstream of the stator and at the exit edge of the stator is calculated by equa- 
tion (2- 10). 

The deviation of pressure rise Apg, through the stator from its time- 
mean value is thus, 



where 

vs 


ir, 


W „ sin 170)0- 

— Exp (ivt) —-Ex 

770) I 17 0) 0“ 


p(il70)0') 


(2-19) 


^PS 

p 



Substituting equations (2-14), (2-15) and (2-iq) into the above and neglecting 
the term of gives 


( 2 - 20 ) 


^PS 


1 - Exp (-iro-j^) 


1 - Exp (-ZircTc) 


® j tan Xr ^ ^ C(wr, 0 -r) 1 + iwg C(wg, Og) 

1 - Exp (-Tr(rg) 


S(ws, <rs) Exp j^-iv(d + is cos Xs)/Va^oj 


rsinjr^^ - Exp .g) I 

1 - Exp L ^ ® J 

S 1 + i.^ "zo - Vu3o tan Xr j ^ 


Exp |^-iv{d +— cos Xs)/Vao Exp (i^t) 
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The dimensionless pressure change is then 


( 2 - 21 ) 


Apg .1 1 - Exp (-irorjE^) 1 - Exp(-2TrtT3) r -i 1 _ Exp (-ttctc ) \ V„ 

7 1 ■ • Fv “ \ , , . 'pr, i S(wc, (Tc) Exp -iv(d + i q cos >l-)/V + tan X-q \ : — . > ^ 

* £-u 2 ^1+ C(wj^, <Tj^) 1 + iwg C(ug, cTg) o S 1^ S S aoj S i + x^>g C(ug, <rg)/ u 


+ tan 


1 - Exp (-trcr } 


1 + iwp '^r) 


r sin TT ug 0'S "I 


- Vu,J 


TT <og trg 


Exp 


r 

|^-iv(d + — . 


Xg)/V3 




Exp (ivt) 


The total ixnsteady pressure difference between the delivery and the 
suction ports of the pump is obtained by summing up hc» ^R> and hg: 


h = 


Pd - Ps /Pd - Ps 




/Pd - Ps \ 
^ £. ..2 yo 


£. 2 
2 ^ 


he + hj^ 


+ hr 


( 2 - 22 ) 


In the case of f = 0, the above becomes 

f 

hqs =: - 2 I tan 1 - Exp (-tto-^)] [ 1 - Exp (-'iro'g) — 

+ tan Xg [ 1 - Exp (-Trcrg)] — 

and agrees, of course, with the pressure change when the axial velocity 
increases steadily by AVg^, assuming the flow through a cascade of flat plates 
without losses. 

Using the result obtained, a sample calculation is carried on for a typical 
single-stage axial flow pump. The basis of the calculation and the result are 
both illustrated in Figure 25. As seen from this result, hg has a phase delay 
that increases nearly proportional to the reduced frequency. The reason for 
this delay is that a part of this pressure fluctuation is caused by the sinusoidal 
gust oscillation produced by the rotor, and the gust requires time to reach the 
stator where the disturbance velocity is converted into pressure. 

It is also worth mentioning that the pressure change due to conduit effect, 
hQ, increases proportional to the frequency, f, and thus becomes predomi- 
nantly larger than the pressure change through the rotor and stator, hg, 

in the higher frequency range. When the frequency exceeds a certain limit, 
the compressibility effect of the working liquid and the piping system becomes 
apparent and the above analysis can no longer be applied. 
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Simplified Dynamic Characteristics of Turbopumps 


Although there are a wide variety of principles of turbopump construction, 
it is fundamentally possible to calculate the dynamic characteristics of indi- 
vidual turbopumps by applying the method developed in the previous paragraph. 
However, the more complicated the pump, the more difficult, and perhaps 
even impractical, is the analytical calculation. In this paragraph, a more 
simple and practical calculation of the dynamic characteristics will be devel- 
oped for the case when approximate dynamic characteristics are sufficient. 

In Figure 26, the shape of turbopumps is different in each case. In all 
cases, however, all of the energy input to the working liquid is produced in 
the rotor. A part of this input is converted into pressure in the rotor, and 
the rest is converted into pressure in the following stator or diffusor. 

For the calculation of simplified dynamic characteristics, the following 
assumptions are made; 

1. The angular momentum of inlet flow into the rotor is steady. 

2. The increase of velocity head at the outlet edge of the rotor is con- 
verted into pressure rise in the following stator or diffusor without any time 
delay. In other words, the stator or diffusor is very narrow and follows very 
close to the outlet of the rotor. 

Under the above assumptions, the change of pressure rise through the rotor, 
App^, is given by equation (2-12) 


^PR 

P 


= - Wu2„AVa 


tan Xr [ 1 - Exp (-irorj^)] 

1 + icoj^ C(wp^, (Tj^) 


Exp (ivt) 


where o-r^ cjR are the stagger angle, solidity, and reduced frequency of 

the rotor cascade, respectively. In the case of the radial and the mixed 
flow impeller, the cascade on the representative stream surface is no longer 
a pure linear cascade, as it is in the case of the axial flow pump. Therefore, 
the result obtained for the axial flow stage does not hold for pumps of other 
types in the exact meaning. 

It is, however, quite reasonable to replace this non-linear cascade 
through an equivalent linear cascade. For example, an equivalent linear 
cascade of a circular cascade (radial flow impeller) can be obtained by trans- 
forming the circular cascade conformally into a linear cascade using a proper 
mapping function. For radial and mixed flow pvimps, therefore, the cascade 
parameters, o*^, and taken not from the original cascade but from 

the equivalent linear cascade. In the same way, axial velocity V , means 
meridian velocity in this case. 

From the previous assumption, the change of pressure rise through the 
stator is equal to the increase of dynamic pressure at the exit of the rotor; 
thus. 
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tan \p [ 1 - Exp (-TTO-p)] _ ,. 

s - V„ AV_ fLi BA! Exp (ivt) 

Uo n 3- jr \ / 


1 + icop_ C(wj^, o-j^) 


The total change of pressure becomes 


^Pr + ^Ps 


= - u • AV - 


tan \p [ 1 - Exp (-tto-d)] 


1 + (Tj^) 


Exp (ivt) 


or by making it dimensionless, 


. _ Ap„ + Apo tan \„ [ 1 - Exp (-ircTp)] AQ 

^R+S - P = - 2 — — 9= ^ — 4, Exp (ivt) (2-23) 

£ 2 1 + 1“R C(03j^, 0-j^) Qq 

2 


.R 


where 


VSR _ cos \p cos \r £_ 
21:11130 Nj^ (j) 4* ^ 


In the case where the solidity of the rotor cascade, <Tj^, is comparatively 
large, namely about o-j^ = 1.0, the above can be simplified further 


r(l/2 + iwj^) AQ 

hR+s = - 2 tan +icop) ^ ^^P 


(2-24) 


r(l/2 + ito) 

The numerical value of is listed in Table I. 

\fn r(l + ico) 

In Figure 27, ^R + g/2 tan 4* i® illustrated in the form of a polar 

diagram. In the case of f = cOp = 0, equation (2-24) becomes 


(% + s)q 3 = - 2 tan — 4) Exp (ivt) 


and agrees with the result when the flow rate changes quasi- statically. The 
ratio, hR + S/(^R +S^qs “ P'(l/2 + icjj^)/'v/iT r(l+io)j^), is also plotted in Figure 28 
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as the Bode diagram, which shows the frequency response by using gain (db) 
and phase shift. From this diagram two limiting reduced frequencies, 

= 0.3 and = 30, are introduced. At a reduced frequency lower than 

(cor)l, the characteristics of a turbopump agree practically with that of steady- 
state operation; that is, the relation between flow rate and pressure rise can 
be expressed by the usual steady- state characteristics curve. On the other 
hand, at a reduced frequency higher than fluctuation of pressure rise 

becomes almost negligible even though the flow rate fluctuates periodically. At 
such a high frequency, the pressure difference due to the oscillating flow in a 
conduit, h^, becomes extremely large. 


Since cor = 


cos X.D f 


Nr«|> 


n 


the corresponding absolute limiting frequencies. 


fj^ and fpj, become 

fL = 0-3 n 

cos \r 

Nd 4^ 

fxT = 30 n = 100 fj 

^ cos \r ^ 


(2-25) 


In the previous analysis, the dynamic characteristics of turbopumps are 
described in the form of frequency responses. It is fundamentally possible to 
calculate the response of pressure rise for an arbitrary change of flow rate 
based on the frequency response. 

The time constant of turbopiimps can be calculated roughly by approximat- 
ing the response of Figure 28 to that of the first order lag element (transfer 
function is K/(Ts + 1)) and by assuming (cor)l =0*3 represents the break point. 
The result is, 


Time constant. 


T = 


0.3 


cos 

2 it Nr (j) n 


(2-26) 


The above shows only the pressure change through a single-stage turbopump. 
In the case of a multi-stage turbopump, the overall dynamic characteristics 
are calculated by summing up the characteristics of each stage. 

Stover all ^ + S 

stage 
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SECTION III. EXPERIMENTAL STUDY OF 
DYNAMIC CHARACTERISTICS OF A CENTRIFUGAL PUMP 


The dynamic characteristics of turbopumps were studied experimentally 
in this Section* The experiment was conducted specifically using a centrifugal 
pump that was available from the commercial market for common industrial 
use. The fluctuating flow was generated mechanically by a reciprocating 
piston attached to the delivery pipe line of the p\imp. 

The instantaneous flow rate and pressures were measured and recorded 
on oscillograph and magnetic tape, which enabled both visual and computer 
evaluation to be made of the measurements. The test result was compared 
with the theoretical prediction that was calculated based on the analysis of the 
previous Section. 

This study was limited to the incompressible case. The experiment was, 
therefore, carried for the frequency range in which the compressibility effect 
is negligible or of secondary importance. Since the occurence of cavitation 
seriously deteriorates the assumption of incompressibility, cavitation in the 
test pump was suppressed as far as possible. 

Experimental Arrangement 

Test Equipment - The test equipment consisted of three main systems; a pump- 
ing system, a pulsating system, and a measurement/control system. Figure 
29 shows the three-dimensional arrangement of the main hardware. 

Pumping system - A centrifugal pump for common industrial use was selected 
for the test pump because of its simple, robust construction and popxilar use. 
Table II describes the details of the test p\imp. 

The pump was driven by a 30 HP induction motor through a continuous 
speed change gear. The moment of inertia of this driving system was so large 
that the fluctuation of input torque to the pump during unsteady operation caused 
no detectable change of rotational speed, so long as the frequency of fluctuation 
exceeded 1 Hz. 

The working liquid (city water) was fed to the pump from a head tank, which 
was located 10 m (33 ft) above the pximp level, contained 26 m^ (7, 000 gallons) 
of water, and supplied 11 N/cm^ (16 psi) boost pressure to the pump inlet. This 
boost pressure corresponded to Net Positive Suction Head of 21 N/cm^ (30 psi) 
and allowed the pump to operate cavitation-free even under the severe test con- 
dition of large pressure fluctuation. 

Suction and delivery acciimulator s were attached to the suction and delivery 
pipe line, respectively. The cushioning action of their air reservoirs served 
to absorb the flow fluctuation in the pipe line. Therefore, the flow fluctuation 
generated by the pxxlsating system could be restricted between the two accumula- 
tors, avoiding undesirable pressure surge in the long suction and delivery pipe 
line s . 
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The flow rate was regulated by adjusting a control valve. In the case of 
fluctuating flow rate, only the mean flow rate, Q^, could be controlled by this 
valve. The distance between the pump and the control valve was designed to 
be as short as practically possible in order to obtain highest possible reson- 
ance frequency. This resonance frequency determined the upper limit of test 
frequency range from which any reasonable incompressible result can be 
reduced. Figure 30 illustrates the detail of test arrangement around the pump. 

Pulsating system - Fluctuating flow through the pump was generated by a 
3 1 /4-inch diameter reciprocating piston attached to the delivery pipe between 
pump outlet and control valve, as shown in Figure 30. 

The instantaneous flow rate, Q, increases as the piston pulls out, while 
Q decreases as the piston pushes in. The fluctuation of flow rate, Q - Qq, is 
fundamentally proportional to the piston velocity. The piston was connected 
to a hydraulic actuator, whose motion was electronically controlled to trace an 
exact sine wave with desired amplitude and frequency. Therefore, the fluctuat- 
ing flow, Q - Qq, is also sinusoidal unless some other deteriorating effects,^ 
such as compressibility, lack of hydraulic power etc., distort the wave form. 

A 20 hp hydraulic pump was used for feeding the actixator with 1,700 N/cm^ 

(2, 500 psi) constant pressure hydraulic oil. 

Measurement and control system - Figure 31 illustrates schematically the 
main measurement and control system. Suction pressure, Pg, delivery pres- 
sure, P(j, and differential pressure, p^ - pg, are measured by variable- 
reluctance type pressure pickups. The deflection diaphram of these pickups 
is interchangeable, so that a diaphram of proper rating is built in for each 
test pressure range. The natural frequency of these pickups ranges between 
15 and 25 kHz, according to their diaphram thickness, and is obviously far 
above the frequency range which is of main interest in this study. 

Flow rate is measured simutaneously by two different kinds of flowmeters, 
a turbine flowmeter and a magnetohydro dynamic (MHD) flowmeter. A 2-inch 
diameter turbine flowmeter was in installed on the delivery line between the 
pulsator and the control valve, and aimed to meter accurate flow rate imder 
steady- state operation. It could not, however, respond quickly enough to 
indicate the instantaneous flow rate under unsteady operation because of its 
relatively slow Digital to Analog converter, which indicated rather the mean 
flow rate, Qq, in a higher frequency range, f > 50 Hz. 

A homemade MHD flowmeter was installed on the delivery line between 
pump outlet and pulsator. This flowmeter consists of an acrylic pipe, a 
magnet of DC exitation and two gold-plated electrodes. The output potential 
from the electrodes was directly proportional to the instantaneous flow rate, 

Q, provided that the flux density across the pipe remains constant. The main 
advantage of this flowmeter was its very quick response, which was of special 
importance in this study. One problem associated with this method was the 
drift of output potential due to the asymmetry of the electro-chemical potential 
around both electrodes. This drift made it difficult to measure the flow rate 
continuously over a long time span with sufficient accuracy. In this experiment, 
however, the MHD flowmeter was operated continuously for about 20 seconds 
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with practically negligible drift. The rate of flow fluctuation, AQ/Qq = (Qmax 
- Qmin)/^Qo> estimated directly by observing the output of MHD flowmeter 

displayed on the screen of the syncroscope standing by the test pump. 

To check the magnitude of the mechanical vibration of the piping system, 
a small piezoelectric accelerometer was cemented on the surface of the delivery 
pipe between the pulsator and the control valve. This measurement was used 
to detect the mechanical resonance of the piping system as well as to detect the 
hyraulic resonance of liquid contained in the piping system. 

The displacement of the actuator plTxnger was also measured by using a 
potentiometer that was built in the actuator. This measurement checked the 
actual motion of the pulsator piston which was supposed to be sinusoidal. 

Seven electrical signals, each representing one of the above measurements, 
were fed to the main amplifier unit, and were recorded both on an oscillograph 
and a magnetic tape; the former for the conventional visual observation, the 
latter for evaluation and processing by digital computer. Necessary calibration 
was conducted to determine the overall conversion factor of each measurement. 

The rotational speed of the test pump was set to a desired value by adjust- 
ing the gear ratio of Varidrive transmission. The rotational speed itself is 
measiired by a digital counter that counts directly the number of revolutions 
per second. This counter also determines the exact flow rate by coxinting the 
pulse signals from the turbine flowmeter. 

A DC power source with automatic current control supplied a constant 
current of 20 A to the magnet of the MHD flowmeter. 

The control system of the actuator consisted of a servomatic analyzer and 
a servo amplifier. The servomatic analyzer generated a sine wave signal of 
the desired intensity and frequency, while the servo amplifier intensified this 
signal and fed the control current to the regulator valve of the actuator. The 
feedback control minimizes any deviation of the actuator motion from the pre- 
determined one. The attenuation of the initial signal generated by the servo- 
matic analyzer was adjusted so that the rate of flow fluctuation, AQ/Qq, was 
a desired value. 

Test Program - An experiment was conducted in accordance with the test pro- 
gram described in Table III. Series 1 and 2 include preliminary tests that 
were intended to determine the nature and magnitude of the pressure fluctuation 
inherent to this test arrangement without applying a forced fluctuation with the 
pulsator . 

The main test was performed at three different rotational speeds of the 
test pump, n = 60, 30 and 15 rps (series 3 through 14). At each rotational 
speed, four series of tests were conducted. First, a steady- state characteris- 
tics curve was determined by varying the flow rate Q from zero to maximum 
without applying any forced fluctuation. (The dynamic characteristics can be 
discussed only in connection with the steady-state characteristics.) 
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Pulsating tests were performed at two different mean flow rates, Q^, 
corresponding to the design flow rate (maximum efficiency point) and to half 
of the design flow rate, respectively. The rate of flow fluctuation, ZSQ/Q^, 
was adjusted to 0. 1 so long as the pulsating system could realize this condition. 
Otherwise, the fluctuation jr ate was maintained at the maximum available value. 
In most cases, the maximum fluctuation rate decreased rapidly when the fluc- 
tuation frequency exceeded 30 Hz. The maximum frequency available from this 
pulsating system was slightly over 100 Hz, but only with a very small fluctua- 
tion rate. 

A linearity test was added to determine the effect of the fluctuation rate, 
AQ/Qq, on the whole phenomena, varying the fluctuation rate from 0.02 to the 
maximum available value at relatively low frequencies, 3 and 10 Hz. As the 
last test, series 15 determined the pressure difference between the pump suc- 
tion and the delivery due to the oscillating flow generated by the pulsator with- 
out running the test pump. This pressure difference, Apc, was essentially 
indifferent to the pumping action itself and merely resulted from the conduit 
effect of the pump casing. This test provided information that was necessary 
to divide the change of pressure difference between the pump delivery and the 
suction, Ap = (p^ - pg) - (p^ - Pg)o> into two essentially different groups, that 
is, into the pressure difference due to the conduit effect of the pump, ApQ, 
and the change of pressure rise through the pump, Apj^^g. 

An oscillograph recording was made for all test series, but a tape record- 
ing was made only for the test series that required computer processing. 

Test Results 

Steady-State Characteristics - From the test series 3, 7 and 11, steady- state 
characteristics of the test pump at n = 60, 30 and 15 rps were obtained. The 
pressure rise through pump, p^ - p^, and the flow rate, Q, were made non- 
dimensional using the following relations; 


pressure rise coefficient. 


4^ = (Pd - Ps>/ 2 




^ (3-1) 


flow coefficient, 4> = Q/tt d 2 b 2 U 2 J 

where d 2 , b 2 , U 2 are outer diameter, exit width, peripheral velocity of pump 
impeller , respectively. 

Figure 32 shows the nondimensional steady- state characteristics of the 
test pump. The design point of this pump was located at <j) = 0. 1137 (=<j>j)) 
where the efficiency of the pump reached its maximum value. As previously 
mentioned, pulsating tests were performed both at cf) = <})£) and cj) = 1/2 

The slope of the characteristics curve is read from Figure 32 as; 
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at (j) = 


d(j) 


5.94 


= ~ 3.44 



The change of pressure rise, ^Pr 4 -S> corresponding to a change of flow 
rate, Q - Qq, can be determined from the above slope when the flow fluctua- 
tion is quasi- static. 

Pulsating Test 

Data and Their Reduction - Pulsating tests were conducted in test series 4, 5, 

6 at n = 60 rps, in series 8, 9, 10 at n = 30 rps, and in series 12, 13, 14 at 
n = 1 5 rps. 

Figure 33 shows a typical oscillograph recording as a sample that illus- 
trates the amplitude- time histories of seven measurements, from bottom to 
top, suction pressure, Pg, delivery pressure, pd, differential pressure, pd 
- Pg, flow rate Q by MHD flowmeter, displacement of the actuator plunger, 
acceleration of the delivery pipe, and flow rate Q by turbine flowmeter. This 
test was a part of series 4 and was conducted imder the condition n = 60 rps, 

Qo = 200 gpm (cj) = cj)j 3 ), f = 10 Hz, AQ/Q^ 0. 1. 

The time passes from right to left on the oscillograph record and each 
abscissa represents a particular time, which can be interpreted from the time 
mark recorded at the topmost edge. This time mark is also simultaneously 
recorded on the magnetic tape to synchronize both recordings. 

Despite the effort to obtain a flow fluctuation as perfectly sinusoidal as 
possible, the one indicated by MHD flowmeter showed in some cases an appre- 
ciable deviation from the initial intention. This deviation resulted partly from 
the compressibility of the water column, and partly from the nonharmonic 
motion of the pulsator piston, which had a tendency to be quicker during pull-out 
motion than during push-in motion, especially in a higher frequency range. The 
pressure fluctuation was also not a simple harmonic but contained many higher 
harmonics that represented either multiple harmonics of basic pressure fluc- 
tuation or standing vibrations of the water column enclosed in the piping system. 
The latter could not be eliminated even with the perfect excitation of fluctuating 
flow. 

Because of the presence of these harmonics, it became practically impos- 
sible to evaluate the correct amplitude and phase correlation between flow rate 
and pressure rise through the visual observation of oscillograph records. 
Nevertheless, oscillograph recording is the best means for visualization and 
understanding of the whole pheonomena involved. 

Test results were recorded also on the multi-channel magnetic tape for 
computer processing. For each test condition, the slice time of recorded tape 
was determined through the visual observation of corresponding oscillograph 
records, so that each slice could contain representative measurements for 
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4 seconds, excluding any measurement under transient or irregular circum- 
stances. First, each slice of the tape was fed to an Analog-to-Digital conver- 
ter, in which the input data on the tape were digitized at the sample rate of 
500 per second to furnish the digital input data for the computer. An IBM 7094 
computor was used for the processing, operating with the Random Vibration 
Analysis Program (RAVAN, Reference 15). The program was developed by the 
Computation Laboratory of George C, Marshall Space Flight Center, NASA, to 
perform various statistical analyses of random processes normally associated 
with flight and captive tests and is capable of computing various statistical 
functions, such as auto- correlation fimctions, power spectral density distri- 
bution functions, root-mean-square (rms) amplitude distribution functions, 
probability density distribution function of single input functions; cross-correla- 
tion functions for positive and negative lag time, cross power spectral density 
distribution functions, cross phase correlation fimctions, and coherence distri- 
bution fxmction of two independent input functions (Reference 16). These com- 
puted functions are available in the form of graphs prepared by a Stromberg- 
Carlson 4020 automatic plotter. 

Figure 34-1 through 34-18 show a part of the plotter outputs for the test 
condition that is equal to that of Figure 33. Flow rate Q by MHD flowmeter, 
change of pressure difference, Ap, suction pressure, Pg, and delivery pres- 
sure, pd, were analyzed as a single input function, and the result is shown in 
Figure 34-1, -2, -3, -4, -8, -9, -12, -13, -14 and - 18. On the other hand, cross 
correlation analyses were conducted for two couples of input functions, (Ap and 
Q) and (pg and p^j), as shown in Figure 34-5, - 6, - 7, -10, -11, -15, -16, and - 17. 

From these resiolts the amplitude and phase correlation of flow and pres- 
sure fluctuation can be determined far more exactly than can be done by the 
visual observation of oscillograph recordings. Besides the accuracy, the 
computed coherence allows the estimation of the reliability of that accuracy. 

Linearity Test - Figure 3 5 describes the amplitude and the phase correlation 
of flow and pressure fluctuation. The flow fluctuation, Q - Qq, is expressed 
by the real part of a rotating vector, AQ Exp (ivt), and the phase of pressure 
fluctuation is correlated to that of flow fluctuation. The change of pressure 
difference in the quasi- static operation condition, Ap^g, can be determined 
from the slope of the steady- state characteristics curve and has 180° phase 
delay to that of flow fluctuation. The fluctuating pressures , Ap, Ap^g, App^ ^ g 
etc. , were made nondimensional by dividing through the corresponding dynamic 
pressure, P^ 2 ^/ 2 y thus giving h, h^g, hp^ ^ S 

Linearity tests conducted in series 6, 10 and 14 were to determine the 
effect of the magnitude of flow fluctuation, AQ/Q^, on the response of pressure 
difference. The result shown in Figure 36-1 (n = 60 rps), 36-2 (n = 30 rps) 
and 36-3 (n = 15 rps) indicates that the nondimensional change of pressure 
difference, h, has the amplitude proportional to the rate of flow fluctuation, 
AQ/Qq, and the constant phase shift over the tested range of AQ/Q^. Since 
the hqg is also proportional to AQ/Q^, the above result means that the ampli- 
tude and phase correlation between h and hqg remains constant, indifferent to 
the magnitude of flow fluctuation. 
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These linearity tests were conducted at relatively low frequencies, at 
f = 3 and 10 Hz. At the higher frequency range where the compressibility of 
the water column becomes significant, good linear relationship cannot be 
expected as in the incompressible case. 

Dynamic Characteristics Test - Figures 37 through 42 show the amplitude 
ratio and phase shift of h/h^g as the function of fluctuation frequency, f, at 
six different test conditions. Because the pressure difference due to conduit 
effect, ApQ, increases proportional to the frequency, the amplitude of h/hqg 
increases also with the increase of frequency. The phase shift varies from 
zero at f = 0 to around 90*" at the frequencies where the ApQ becomes pre- 
dominantly larger than ApR 

The effect of compressibility becomes appreciable when the frequency 
exceeds 30 Hz. The intensity of pressure fluctuation increases rapidly as 
the frequency becomes close to the resonance frequency of the pumping 
system, around 120 Hz. Figure 43 illustrates the change of pressure differ- 
ence, Ap, corresponding to the flow fluctuation of AQ = 10 gpm using a polar 
diagram. The pressure difference due to conduit effect, ApQ, was determined 
from the pulsating test of series 15 as follows; 

amplitude of ApQ (psi) = 0.0277 AQ(gpm) f(Hz), 

phase shift of Apc =96® delay to flow fluctuation. 

The conduit effect of the pump casing is equivalent to that of a 38.4-inch long 
pipe with a 2-inch inner diameter. 

Subtracting Ap^ from Ap vectorially, the change of pressure rise through 
the pump, App^^g, can be determined. Then, the amplitude and phase correla- 
tion between stnd hqg can be established. 

The amplitude and phase correlation of ^R + g/hqg are plotted in Figures 
44-1 and 44-2 against the reduced frequency of the rotor cascade, 

_ cos ^ 

“R “ Nr 4) n 


where =55® and Np^ = 6. 

The data shown in these figures were confined to those that seemed to have no 
compressibility effect. The test results were relatively scattered in spite of 
very careful reduction of data by computer. This was mainly because App^g 
cannot be determined directly but indirectly by subtracting Ap^ from Ap 
vectorially. Especially in higher frequency range, App must be determined 
as the difference of two vectors, Ap and Ape, that are predominantly larger 
than App^ ^g itself, and the accuracy of obtained App ^g inevitably deteriorates. 

To compare the test result with the theory, a theoretical curve was drawn 
which was calculated from equation (2-24) of Section II assuming a simplified 
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model of turbopumps. The test results showed the tendency for the amplitude 
ratio of h|^ +s/^qs 'to decrease with the increase of the reduced frequency at 
a slightly faster rate than the simplified model theory predicted. Regarding 
the phase shift of hj^ + s/^qs» the experimental shift was definitely larger than 
that of the simplified model theory. This suggests that the conversion of 
kinetic energy at the exit of the impeller into the pressure rise in the following 
volute chamber requires a certain time delay that is neglected in the simplified 
model theory. 

It is, of course, possible to calculate the frequency response of the test 
pump more accurately by considering the mechanism of pressure recovery in 
the volute chamber, thus reducing the discrepancy between theory and experi- 
ment. Nevertheless, the simplicity and the generality of the simplified model 
theory seems to be worth even more than its reduced accuracy. 

The reduced frequency, a-t which the dynamic characteristics 

begin to separate from the quasi-static ones, is read from the test results as 


= 0.1 


comparing with =0.3 obtained from the theory. 

The time constant of the response of pressure rise to a step-like change 
of flow rate was calculated by approximating the obtained frequency response 
to that of a firstorder lag element and by assuming =0.1 represents 

the break point. The result is 


Time cons ant, T 


1 cos 
0.1 Ztt <\> n 


(3-2) 


Regarding the relation between mean flow rate, Qq, and corresponding 
mean pressure rise, (p^ - Ps)o» during the fluctuating operation, it is found 
that the relation is practically identical to that of the steady- state character- 
istics curve as far as the rate of flow fluctuation, AQ/Q^, does not exceed 
the test range. 

Fl uctuation of Suction and Delivery Pressure - According to the force flow 
fluctuation, suction pressure, p^, and deliver y pr e s sure , p^, fluctuate with 
the same frequency. The amplitude ratio of pj and Pg is plotted against fre- 
quency in Figure 45. In the lower frequency range, the amplitude of pj is 
much larger than the amplitude of Pg, because the pressure fluctuation at the 
suction port is suppressed by the presence of the suction accumulator. At 
higher frequencies, however, both amplitudes become more or less equivalent. 
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The phase correlation of and pg is shown in Figure 46. Since the whole 
test frequency range is under the resonance frequency, both pressure fluctua- 
tions have comparatively less phase difference. 

In this study cavitation in the test pump was prevented by supplying a large 
boost pressure at the suction. In the case in which the boost pressure is not 
large enough, the fluctuating suction pressure often causes periodic cavitation 
which introduces strong nonlinearity to the system. Under the cavitating con- 
dition, it was observed that the frequency of the predominant pressure fluctua- 
tion was not the frequency of forced fluctuation but of subharmonics, f/2, f/3, 
f/4 etc. 


CONCLUSIONS 

The response of pressure rise of turbopumps to periodically fluctuating 
discharge, termed dynamic characteristics, was studied analytically and 
experimentally. 

In Section I, the theoretical basis for the calculation of dynamic charac- 
teristics was developed. The incompressible, inviscid flow solution around 
a cascade of nonstaggered flat plates for a periodically fluctuating inlet flow 
was determined. Solutions were obtained for two kinds of inlet unsteadiness, 
namely for translatory oscillation and for sinusoidal gust oscillation. 

The possibility of applying the results of the nonstaggered cascade to the 
cascade of arbitrary stagger angle was sought. This analysis proved that the 
obtained solution was applicable also to arbitrary cascade with sufficient 
accuracy. 

In Section II, the analytical method for calculating the dynamic character- 
istics of turbopumps was developed, applying the results obtained in Section I. 

A single-stage axial flow pump was chosen as the sample case and its exact 
dynamic characteristics were determined based on the specific data of the 
rotor and stator cascade. 

For the case when even approximate dynamic characteristics are sufficient, a 
simple method was also derived that enables the calculation to be made without 
entering into the details of the pump design. 

In Section III, the dynamic characteristics of a centrifugal pump were 
determined experimentally. Pressure at suction and delivery, flow rate, etc. 
were measured and the correlation between pressure and flow fluctuation was 
analyzed. 

The dynamic characteristics obtained from the experiment were close to 
those predicted by the simplified model theory of Section II. The limiting 
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frequency under which the dynamic characteristics closely agreed with the 
quasi- static ones, was determined experimentally as 


f = 0. 1 


cos 


n 


In the entire study, interest was confined to the incompressible, non- 
cavitating case. When cavitation occurs in a turbopump, the phenomena change 
drastically from those described in this work. 
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TABLE I. value OF V (1/2 + icj)/ nJtt F (1 + iw) 


r(i + ico) 

\fn~ r( 1 + ico) 


R + il = 


CO 

R 

— 

I 

P 

e’ 

0 

1 

0 

1 

0 

0. 01 

1. 000 

-0. 014 

1. 000 

- 0. 79 

0. 03 

0. 998 

-0. 041 

0. 999 

- 2. 38 

0. 1 

0. 975 

-0. 134 

0. 984 

- 7. 81 

0. 2 

0. 910 

-0. 242 

0. 941 

-14. 89 

0. 3 

0. 826 

-0. 314 

0. 884 

-20. 80 

0. 5 

0. 669 

-0. 372 

0. 764 

-29. 04 

0. 7 

0. 554 

-0. 371 

0. 666 

-33. 75 

1.0 

0.447 

-0. 342 

0.563 

-37. 44 

1. 5 

0. 352 

-0. 297 

0.461 

-40. 12 

2 

0. 299 

-0. 264 

0. 402 

-41. 38 

3 

0. 240 

-0. 221 

0. 326 

-42. 60 

5 

0. 183 

-0. 174 

0. 252 

-43. 57 

10 

0. 127 

-0. 124 

0. 178 

-44. 28 

30 

0. 074 

-0. 073 

0.103 

-44. 77 

100 

0. 042 

-0. 042 

0. 056 

-45. 00 

00 

0 

0 

0 
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TABLE II. SPECIFICATION OF TEST PUMP 

Model / Manufacturer 

1 CCL / Buffalo Pumps, N.Y. 

Suction / Discharge Diameter 

2 inch / ^ “ inch 

Impeller Type 

Radial Impeller with Backward Vanes 

Outer Diameter 

6 , 6'32 inch 

Exit Width 

0.260 inch 

Number of Vanes 

6 

Exit Angle of Vanes 


Material 

Bronze 

Diffusor Type 

Vaneless Volute Chamber 

Casing Material 

Cast Iron 

Normal Operating Condition 


Rotational Speed 

3,500 rpm (38.33 rps) 

Flow Rate 

120 gpm 

Head Rise 

130 ft 

Power Required 

10 HP 



TABLE ni 


TEST PROGRAM 


Series , • ^ ^ „ 

actuating 

test pump system 

valve 

settings 

measurement 

NO. 

of test 

n Qq hyd . f 

(rps) (gpm) ^Q/Qq pump (Hz) 

main items 

record on 
O” tape 
gram 

1 

, , _ AV: open & 

preliminary test to check pressure u q ^ run ^ closed 

fluctuation caused by hyd. pump (stop) (2500. v CV: varied 

_ . _ . -n gj l 

Ps’ Pd’ Pd-Ps’ 

^MHn ^ ^turb. 

X 

2 

preliminary test to check pressure ^ , 0 to 

fluctuation at min- disturbance max. 

0 stop 

DS: closed 
0 DV:open 

CV; varied 

ditto 

X 


3 

0 to 

steady-state characteristics 60 max. 

0 stop 

0 regular 

ditto 

X 


4 

pulsating test at c|) = 200 

0.1 or 
, run 

less 

1 to 

regular 
100 ^ 

ditto 

X 

X 

5 

pulsating test at (|) = •” 100 

0.1 or 
, run 

less 

1 to , 

100 

ditto 1 

X 

X 

6 

linearity test 60 200 

0.02 tc 

run 

max. 

3 & 1 C regular 

ditto 

X 

X 

7 

steady-state characteristics 30 ^ 

^ max. 

0 stop 

0 regular 

ditto 

X 


8 

pulsating test at (|> = 30 100 

0.1 or 
less 

run 

1 to 

100 regular 

ditto 

X 

X 

9 

pulsating test at = •“ 4 ^^ 30 50 

0.1 or 
less 

run 

1 to 

100 

ditto 

X 

X 

10 

linearity test 30 100 

0.02 tc 
max. 

run 

3 & 10 

regular 

ditto 

X 

X 

11 

steady-state characteristics 

0 to 
max. 

0 

stop 

0 

regular 

ditto 

X 


12 

pulsating test at 4* = 

15 

50 

0.1 or 
less 

run 

1 to 
100 

regular 

ditto 

X 

X 

13 

pulsating test at ^ 

15 

25 

0.1 or 
less 

run 

1 to 
100 

regular 

ditto 

X 

X 

14 

linearity test 

15 

50 

0.02 tc 
max. 

run 

3 & 10 

regular 

ditto 

X 

X 

15 

pressure fluctuation caused by 
actuator motion only 

0 

0 

CD 

run 

1 to 
40 

regualr 

ditto 

X 

X 


Regular valve settings; suction stop v. , delivery stop v. (DS) , accumulator v. (AV) open, 

drain v (DV) closed, control v. (CV) adjusted. 
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FIGURE 1. UNSTEADY OPERATION OF A TURBOPUMP; 
FLOW RATE AND PRESSURE FLUCTUATE PERIODICALLY 
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FIGURE 2. UNSTEADY FLOW THROUGH A CASCADE 


FIGURE 3. FLOW CONFIGURATION AROUND A CASCADE OF 
NONSTAGGERED FLAT PLATES 

Airfoil Wake 



FIGURE 4. THREE KINDS OF CIRCULATION FOR 
TRANSLATORY OSCILLATION 
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FIGURE 9. CIRCULATION ENCLOSED IN INTEGRATION PATH 


vortex sheet, 



FIGURE 10. SIMPLIFICATION OF FLOW FIELD DOWNSTREAM 


OF THE CASCADE 
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FIGURE 16. COMPARISON OF CIRCULATION FOR STAGGERED AND 
NONSTAGGERED CASCADES 
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FIGURE 17. EFFECT OF STAGGER ANGLE ON THE CIRCULATION 
OF AN AIRFOIL IN A CASCADE 



FIGURE 18. SINGLE-STAGE AXIAL FLOW PUMP FOR ANALYSIS 
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FIGURE 21. STEADY AND UNSTEADY FLOW FIELD 
AROUND A ROTOR 
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- AVg^ tan Exp (ivt) 

(a) Through-Flow (b) Vertical Flow 


AY 3^ tan Exp (ivt) 
(c) Induced Flow 


Irrotational Flow 

o, 


Rotational Flow 
(^l>does not exist) 





(d) Simplified Induced Flow 
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FIGURE 22. VELOCITY POTENTIAL CORRESPONDING TO EACH 
UNSTEADY FLOW ELEMENT 



Sinusoidal 

Steady Translatory Gust 

Flow Oscillation Oscillation 



FIGURE 23, STEADY AND UNSTEADY FLOW FIELD 
AROUND A STATOR 
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(a) Modification 


FIGURE 24 . VELOCITY POTENTIAL CAUSED BY SINUSOIDAL 

GUST OSCILLATION 
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FIGURE Z5. DYNAMIC CHARACTERISTICS OF A SINGLE-STAGE 

AXIAL FLOW PUMP 
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FIGURE 28. BODE DIAGRAM OF S^^^R + S^qs 
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FIGURE 29. SCHEMATIC VIEW OF THE EXPERIMENTAL ARRANGEMENT 





FIGURE 30. DETAILS OF THE PUMPING AND PULSATING SYSTEM 
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FIGURE 31. SCHEMATIC DIAGRAM OF THE MEASUREMENT 

AND CONTROL SYSTEM 
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FIGURE 33. TYPICAL OSCILLOGRAPH RECORDING OF UNSTEADY- 

STATE OPERATION: n = 60 rps, = 200 gpm, f = 10 Hz, 

AQ/Q =0.1 (SERIES 4) 
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FIGURE 34-1 TO TYPICAL SC 4020 PLOTTER OUTPUTS PROCESSED 
34-18 BY THE RANDOM VIBRATION ANALYSIS PROGRAM 

(RAVAN): n = 60 rps, Qq = gpm, f = 10 Hz, AQ/ Qq s; 0. 1 

(SERIES 4) 
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FIGUBE 34-3.rms AMPLITUDE DISTBIBUTION OF Q VERSUS FREQUENCY 
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FIGURE 34-5. CROSS POWER SPECTRAL DENSITY OF Ap AND Q 
VERSUS FREQUENCY 
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FIGURE 34-6. CROSS PHASE CORRELATION OF Ap AND Q VERSUS FREQUENCY 
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FIGURE 34-8. AUTO- CORRELATION OF Q VERSUS LAG TIME 
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FIGURE 34-9. AUTO- CORRELATION OF Ap VERSUS LAG TIME 
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lag time (unit=2 ms) 


FIGURE 34-11. CROSS-CORRELATION OF Ap AND Q VERSUS NEGATIVE 
LAG TIME 
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FIGURE 34-13. rms AMPLITUDE DISTRIBUTION OF SUCTION PRESSURE, p^, 

VERSUS FREQUENCY 
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FIGURE 34-14. rms AMPLITUDE DISTRIBUTION OF DELIVERY PRESSURE, 

p^, VERSUS FREQUENCY 
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FIGURE 34-15. CROSS POWER SPECTRAL DENSITY OF Pg and VERSUS 
FREQUENCY 
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AMPLITUDE DEVIATION 
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FIGURE 38. RESPONSE OF PRESSURE DIFFERENCE 
OBTAINED FROM PULSATING TEST: n = 60 rps, <p = 

(SERIES 5) 
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FIGURE 39. RESPONSE OF PRESSURE DIFFERENCE 
OBTAINED FROM PULSATING TEST: n = 30 rps, =<#>q 

(SERIES 8) 
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FIGURE 40. RESPONSE OF PRESSURE DIFFERENCE 
OBTAINED FROM PULSATING TEST: n = 30 rps, cp = 4>jy/2 

(SERIES 9) 


101 


PHASE SHIFT OF hA (d^g) 



70 



FIGURE 41. RESPONSE OF PRESSURE DIFFERENCE 
OBTAINED FROM PULSATING TEST: n = 15 rps, 

(SERIES 12) 
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FIGURE 42. RESPONSE OF PRESSURE DIFFERENCE 
OBTAINED FROM PULSATING TEST: n = 15 rps, 0 = 0^/2 

(SERIES 13) 
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FIGURE 43. POLAR DIAGRAM DESCRIBING RESPONSE 
OF PRESSURE DIFFERENCE FOR THE FLOW FLUCTUATION OF 
AQ = 10 gpni (SERIES 4, 5, 8, 9» 12, 13, and 15) 
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FIGURE 45. AMPLITUDE RATIO OF DELIVERY 
AND SUCTION PRESSURE 
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FIGURE 46. PHASE 
AND 
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